Theory of Second Order Optical Processes from A Luttinger 

Liquid 



Jun Zang, Joseph L. Birman 
Physics Department, The City College of The C.U.N.Y 
Convent Avenue at 138th Street, New York, NY 10031 

Zhao-Bin Su 

Institute of Theoretical Physics, Academia Sinica 
Beijing 100080, People's Republic of China 

Abstract 

We develop a theory for the total optical secondary emission from a ID inter- 
acting electron system modelled as a Tomonaga-Luttinger liquid. We separate 
the emission into two parts which may originate in hot luminescence (HL) and 
Raman Scattering (RS) respectively when we neglect the interference effect. 
We hnd a peak around Aw = Vf\q\ in the RS part which does not come from a 
structure factor peak. In general the total emission cannot be separated into 
HL and RS. However at resonance, and taking into account the k dependence 
of the optical matrix element, a part of the RS is proportional to the structure 
factor S{qi — q2,uJi — ^2). 
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Typeset using REVT^ 



Recent advances in fabrication of semiconductors permit the construction of systems, 
such as the quantum wire, with effective geometrical dimensionahty of one (ID). This has 
permitted study of the transport properties of ID systems. The optical properties can 
also be studied, so as to compare and contrast with electronic models (Fermi Liquid or 
Luttinger Liquid) for ID electronic and optical response. It is well known that the Tomonaga- 
Luttinger (TL) model [Q is necessary to describe a pure ID interacting electron system 0. 
However, in an actual semiconductor quantum wire, the scattering and finite size effects 
may serve to damp out the low-energy plasmons and restore the Fermi surface [^. It is not 
yet clear whether the TL liquid or the normal Fermi liquid theory is more appropriate to 
describe properties of semiconductor quantum wires. As a matter of fact, both TL liquid 
and Fermi liquid theory have been used to study the transport and optical properties of ID 
semiconductor electron systems @-§|,PJ31- 

Optical methods are powerful tools for the investigation of the collective and single parti- 
cle excitations. For example, the Fermi edge singularities (FES) 0|1O| in the photoemission 
and soft-x-ray emission and absorption provide information about the existence of the Fermi 
surface and the low energy excitations. The conditions for FES have been studied in both 
normal Fermi liquid theory and TL liquid theory Although Raman Scattering 



(RS) experiments have been carried out in these ID quantum wires |Tl|, there are only a 



few theoretical studies using normal Fermi liquid theory [0 which aim to describe the RS 
spectra of these systems. However, the calculation of the RS efficiency from the collective 
excitations in ID systems needs special caution, especially in the resonance region, as we 
discuss later. The objective of this paper is (1) to give a theory of RS in a TL liquid; (2) 
to reveal that generally the second order RS efficiency is not proportional to the dynami- 
cal structure factor in the resonance region contrary to a widely used approximation; (3) 
to discuss the difference between RS from TL and that from FL. We also remark on the 
possibility of seperating hot luminescence (HL) and Raman scattering (RS). Since TL is an 
exactly solvable model, our calculation may serve as a special prototype for the discussion 
of various aspects of RS, such as threshold behavior. 



Let us first consider the general RS which is due to the A ■ P interaction. Our discussion 
here will follow the standard literature on the RS from collective excitations [^,0. If we 
here neglect the "final state interaction" of the core holes, then in second order perturbation 
theory, the transition rate can be written as: 



dr / dt' / 



dt 



X (0|cL,,(t' + r)c,,_,,(r)cI_^^(0)c,_,,(t)|0)e-^(-+-)(*-*')-^(---)^ (1) 

Here gi, ui (g2, ^^2) are the wave vector and frequency of the incident (scattered) photons, 
(c|,) is the annihilation (creation) operator of the conduction band electron, ed is the eigenen- 
ergy of the core hole, and Mk is the optical matrix element. The quasi-particle approximation 
can be used to simplify the calculation as well as to "pick-up" only the virtual RS process 
contribution. In this approximation, we substitute c^/_qj(t' + r) = e*'''='-«i* c|./_^^(r) and 
Ck-q^{t) = e-*''=-n*Cfc_g(0) m in Eq.Q], and we obtain 

W{qu gs; uJi, UJ2) = 2^ j j 

k,k' ^1 ~ ^d) t^l — [^k'-qi — ^d) 

X /''°'rfr(0|cI_,^(r)c._,,(r)cI,_,^(0)c.,_,,(0)|0)e-^(— )^ (2) 

If we treat Mq as a constant, and assume the RS is near but not at resonance, i.e. incident 
photon frequency Ui ~ Ag but — {^k-qi — q)| >> Cfc-gi, where is the semiconductor 
band gap, then we obtain the celebrated form |T^,[I^ W{qi,q2]UJi,uj2) = (^JJ^j^ ^i S{qi — 
q2,uji — UJ2) and the Random-Phase Approximation (RPA) can be used to calculate the RS 
efficiency from collective excitations. In this case, the RS efficiency is proportional to the 
dynamical structure factor S{qi — q2,uJi — UJ2) with an enhancement prefactor. Although 
this is a very useful result, its applicability in the resonance region is invalid, especially for 
the ID electron gas. 

Since the ID TL model is exactly solvable even including the core hole effect, it would 
be of interest to study RS from the TL liquid. Before including the core hole effect, we 
show that taking into account the /c-dependence of there is a contribution in Eq. |1| 



which is proportional to the structure factor S{qi — q2, uji — U2) for RS from the TL hquid. 
Expand around some point [|I^ and around some point k'j^ {k'j ~ Fermi momentum 
kf). Since in resonance RS, we can choose the appropriate ko and k'^ then ''-II.) ^ 

|Mp 2Re(M'M*) 

aJi-(efc_,j-etj) Vf 



Here f/ is the Fermi velocity. It is easy to see that the 2nd term 



will give a contribution proportional to the structure factor 5(^1 — g2, ^1 — ^2) in Eq[0- The 
relative amplitude of this contribution depends on detailed information about the system. 

Let us now calculate the RS efficiency for the spinless TL liquid with constant matrix 
elements Mk = M. Experimentally the spinless model can be realized by applying a strong 
magnetic field that polarizes electron spins. Generalization to the spin- 1/2 case is straight- 
forward but the calculation is more complicated. In the standard bosonization of the ID 



electron gas, the Hamiltonian of the coupled TL bosons and localized core hole is |16 



Ho = Yl '>^TL\k\hlh + E e<i44 + ^ E Vue'^^dld^h + h.c. (3) 



k X \ k,x 

The electron-radiation interaction can be written as: 



Hi^t = MJ2 ^i{x)d,A^e''^^ + h.c. (4) 

x,q,a 

where ip^ix) is the fermion operator of right moving (a = +) and left moving (/3 = — ) 
electrons, and bl and bk are the operators of the TL boson whose velocity is vtl = 
[{vf + VikY — (V2fc)^]^''^. dx and d'^ are operators of the localized core hole with energy e^. 
Here we only include the forward scattering of the electron by a core hole in Eq. 0, the effects 
of backscattering of electrons will be discussed in the end. is the renormalized interaction 
between the TL boson and the core hole Vk = V3fc(^)^/^[cosh 0^ — sinh 0^]. Here we use the 
notation as in for Vik, and 0^. Vik and are for the forward Fermion scattering in- 
teraction, so in "(yf-ology", Vik = and V2k = 92- Vs^ is the interaction between electron and 
core hole. In the further calculation, we will use the short-range 5-potentials for Vi{x). The 
effect of Vi is a shift of the Fermi velocity, so we will define vj = vj + Vi. For the 5-potentials, 



the effective Fermion- Fermion coupling constant g = sinh {(pk) = |[1/y 1 ~ ^2/^/ ~ !]• To 
calculate the transition rate W{qi,q2;uJi,uj2), we need to calculate an 8-point correlation 
function if we include the core hole effect: 



C(x,x',M,t',t) = (0|^,(x',t')4'(O4'(ti)^l(x',M)^/3(x,0)4(0)4(t)4(x,t)|0) (5) 



Using a unitary transform |17|, the 8-point correlation function decouples to two 4-point 



correlation functions, which can be calculated easily. The transition rate is calculated as: 

W^"^(gi,g2;wi,W2) = const f dye^^^^^^-^^)^ f du f d^;e-*('^2«+-i^) 



X f^y^y^ ^ (6) 

[{u + y — ie) iy — y + ie)]^"'' [{u — y — ie) {v + y + ie)]^"'^ 



With 



\Ax^ - iu- V + 2y -2ief]^"0\Ax^ - iu- V -2y -2ieY]^'-i^ 
f(u,v,y)= I dx — — 7 

Here e is the inverse of the ultraviolet momentum cutoff, and we have renormalized the 
photon frequencies by cUj ^ cijj — (cj — e^). The g's are defined as: g+ = g\j^ = gZ_ = 
[v/^-(5(v/^-y^)]2; g_ = g-^ = g+_ = _ 5 - ^)]2; g^^^ = ^/^^ 

(a 7^ A useful combination is ?7 = g+ + g^. Here 6 = V^/vtl is the phase shift of 



noninteracting electrons scattered by V3 [0 



Our exact results Eq. g-j^] are similar to that in Ref. where Nozieres and Abrahams 
discussed the threshold behavior of RS from noninteracting electron systems. But they did 
not discuss the interference effects. If we neglect the space-dependence and the electron- 
electron interaction effects, we will get the same results as in Ref. |T^. Note however the 
exponents we calculated here are different because of the coupling of the core hole to both 
branches {k and —k) of the electrons. 

Now let us calculate the transition rate W{Q, 001,002) with Q = qi — q2- We drop the 
superscrpt a = (3 = + used in Eq. and write it later. We can obtain the result for 
a = (3 = — hj changing Q —Q in W{Q, 001,002). The calculation for a 7^ /3 is similar 
and we will only discuss the results. It's easy to see there is a divergence in Eq.0|, and this 
divergence can be cured by introducing the core hole lifetime r. In addition to removing 
the unphysical divergence, we introduce the core hole lifetime r, as in Ref. |jl8[ in order to 



be able to separate "hot luminescence" (HL) and resonance Raman scattering (RRS) in the 



total secondary emission given by Eq.p. Following Ref. [|^, f{u, v, y) can be split into two 
parts: f{u,v,y) = T + \\u + v\^!{^i^^,^^ with 

After substituting Eq.|^, the first term in Eq.p] which is proportional to r will give the 
hot luminescence (HL) contribution, and the second term will give the Raman Scattering 
(RS) contribution. Define uj± = uJi — uj2^ Qvtl, for t] > 1, then the HL contribution to the 
transition rate is [O: 



Wf{Q, 001,002) = canst 0(cJi)e(-cj2)(-cu2)^~^e-^('^^-'^2) / dss^—\l - sy+-^ 

Jo 

{2uj2S + oo^y—\2uJi - 2uo2S - uo_y+-^Q{2uo2S + oo^)Q{2uJi - 2uo2S - oo_) (9) 

The numerical calculation of Wp is shown in Fig.(l). The calculation shows that there 
is singular behavior around uoi — 002 = \Qvtl\- Since 002 < 0, so from Q{2oo2S + 00-) and 
0(2ci;i — 2uo2S — uo^) we have WpiQ, 001,002) oc 6(ci;i —002 — \Qvtl\)- The other two peaks are 
at ijOi+ 002 = ±Qvtl- The asymptotic behavior of Wp can be given at these singular peaks 
using the following integral 

/ x"(c - x)%x + a)^ = c'+^+^a^ ^/o n 2Fi{-w, 1 + u,2 + u + v, -c/a) (10) 

JO L [2 + U + V) 

Where 2-^i(— w,l + u,2 + u + v,—c/a) is the Hypergeometric function. Without losing 
generality, we assume here Qvtl > 0. Define c = {ooi — UO2 — Qvtl) I {—'^002) and e = 
{ooi + 002 + Qvtl)/{—2oo2), then e = c + {^z^ ~ !)• Now we discuss various limits. 

(i) at c ~ 0, and < e < Qvtl/ {—002) [as in Fig.(l)] 

Wf-- (-^^2)"""^ r dss^—\c ~ sy'-\s + ey+^^ 
Jo 

~ ^3-h^^-\F,{l~ g^,g.,2g_,-c/e) 

If c < e, then Wp ~ c^^—i- If O e, then Wp ~ Cic^3—9+^'2{9+~i) + Q^^zg--i^9+-i~9~ ^ 
where Ci is constant. If c ~ e, then Wp ~ c^-"'"'^"^. 

(ii) at c ~ 0, and e < 



~ 0(c - |£|)(c - - g.^g^, 2r], 1 - c/\e\) (11) 

If c > then Wp ~ Cic3+ i) + Csc^— If c ~ then - 
(c- 

(iii) at £ ~ QvtlI {—^2) then c ~ 1. Define A = c — 1, then 

JO 

r CiA^f— 2 + CaA^'-i, if A > 0; 

(12) 

i(:7i|A|2f+-2 + C2|A|3+-f— 1, ifA<0 

(iv) In the above cases, we have assumed UJ2 9^ 0. If \uj2\ ^ oj^ and Icui + UJ2 ^Qvtl\ ^ |i^2|, 
then 

Wf ~ (-u;2)''~'(c^-)'-"'(c^i + a;2 + gt^L)"-^"' (13) 

For ?7 < 1, the HL contribution will diverge. We believe this divergence comes from the 
interference effect which we included. The calculation of the RS part contribution for r] <1 
is an extremely tedious task. However since in this case the HL part will diverge and the 
integral of the RS converges well, we will not discuss the calculation in this region. For 
?7 > 1, the integration in Eq.p] can be estimated around the lower limit of the integration. 
The RS part can be calculated: 

Wr{Q,uoi,uo2) = const {ui + U2y'~h^~"e{ui - U2MU-) + 0{e)f-iu-)] 

[<5(^+) + 0(e)/+(o.+)] + --- (14) 

Here f^{ujJ) and /+(ct;+) are smooth functions of uo^ and uj+ respectively. 

We can see that there is a peak around oji — UJ2 = \Qvtl\- The prefactor {ui + uj2y~^ 
clearly shows that this peak has a different origin from the peak of the structure factor. 

In the first order optical process such as photoemission and absorption, the competition 
between electron-electron interaction and electron-hole coupling in the TL is manifest in the 
FES exponent; the condition for the critical point of positive and negative exponents is 77 = 1 
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P]. In secondary optical processes, because of the interference of light from recombination 
of holes and electrons at different sites, the spectrum is more complicated, the competition 
of electron-electron interaction and electron-hole coupling will be manifest in the RS part 
as well as HL part. But the critical point is the same [r] = 1). 

Now we turn to discuss W^^ with a ^ (3, which is related to the processes which describes 
the situation where photon excited electron and the recombination electron are in different 
branches. For 77 > 1, the RS contribution is the same as Eq. [|14|. The HL contribution 
is proportional to T^~^V9^~V9+f << So for y/g~g+ > 1/2, the HL contribution can be 
neglected. If we assume the core hole lifetime r is larger than the characteristic time scale, 
which will cure the divergence when 77 < 1, the HL contribution for ^J'gZgl^ < 1/2 is also 
negligible compared to the a = P case. 

It is tempting to discuss the threshold behavior of RS from a higher-dimension (2D or 3D) 
electron system as Schotte and Schotte did for FES in emission (absorption) |T^. However, 
the possibility of extending our result to 3D depends essentially on the scattering phase shift 
being pure s-wave, which is not true in our case, since we do include interference between 
different core hole states. Technically this prevents carrying through a partial wave expansion 
which is needed for the ID — 3D analogy to be valid from TL liquid absorption/emission to 
threshold optical behavior of the 3D Fermi gas. But we can still use our result to discuss the 
interference effects and interaction effects in Fermi liquid in high dimensions. Following the 
recent arguments in renormalization group theory of Fermi liquid pOf , this is especially true 



in the threshold region of the secondary optical processes. When neglecting the interference 
effects, similar to the absorption (emission) case, the threshold behavior of RS in a high 
dimension (2D or 3D) non-interacting Fermion gas is similar to that of the TL liquid [jl8 



From our above results, one effect of the interference which seems to be true in a high 
dimensional Fermion gas is that Wp can not be factorized into absorption and emission 
any more. Similar effects have been found in RS from polaritons |21[. In Fig.[l], there are 
features around c<Ji — co'2 = \Qvtl\ which suggest Wp looks more like Raman Scattering than 
luminescence. 
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A more striking effect in the TL liquid is that the transition rate is very sensitive to the 
couphng strength of TL bosons to the core hole, and the value of g. We show in Fig. [2] the 
range of V2 and V3 that gives the divergent and convergent HL Wp- We can see from Fig. [2] 
that large g (or V2) and small V3 is needed to have convergent Wp- A interesting result is 
that the Wp diverges for a ID Fermi Liquid {g = 0). So the HL part of contribution will 
dominate the secondary optical processes. The coupling of electron-hole (V3) will increase 
this effect. When the electron-electron interaction is strong enough, the HL part Wp will 
converge and the RS Wr will dominate the total secondary emission. Since the HL part and 
RS part have different peak (singular) structure, the theory does show different secondary 
optical properties between TL and FL. 

Strictly speaking, the HL part and RS part can not be called "Hot Luminescence" and 
" Raman Scattering" separately, since Wp can not be factorized into absorption and emission 
in the usual sense, and there is singular structure around ui — UJ2 = \Qvtl\ in W^. But 
we think of them as having originated from "Hot Luminescence" and "Raman Scattering" 
respectively although the interference effects make them non-separable. 

Our main objective in this Letter has been to give an exact model of secondary optical 
process in the resonance region. For clarity and simplicity, we used a model of a localized 
hole coupling with the TL without backscattering. Also the electron-electron interaction 
is taken as short-range (5-potential). The effects of finite hole mass and unscreened long 
range electron-electron interaction need further study. The most important effect is that we 
omitted the backscattering of the electron by core hole in TL liquid. In the point poten- 



tial, p2| , p3| each correlation function in Eq.(^ can be written as product of a correlation 
function due to forward and backscattering potential G = GjGb- So there are addtional 
contributions to both Wp and Wr. It was shown by a RG calculation in Ref. |Q that a 
TL liquid with repulsive interaction will be decoupled into two separate TL liquids by even 
a weak impurity if the backscattering is taken into consideration. There are recent studies 
of backscattering effects on the FES in photoemission and photoabsorption [^^-[2^. In the 



weak interaction limit, Kane et al. showed that there is a crossover energy between strong 
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and weak backscattering for e* ~ (Vfe/^f j)^/'^^ ^\ At e < e* , has power-law behavior, and 
the exponents due to backscattering are positive. We suggest this result is generally true 



even for strong electron-electron interactions p2|,^. Since the divergence of the HL is due 
to the low energy excitations, the backscattering effect will make the HL converge better. If 
the core hole lifetime is smaller than e*, the HL contribution Wp will be changed mainly due 
to the increase of g±. If the core hole inverse lifetime is larger than the backscattering energy 
interaction of electron and core hole, the backscattering effect on Wp will be less dramatic, 
and we suggest here the main feature of Wp will not be changed. The RS is a virtual process 
and independent of the core hole lifetime. From Eg. ([T^) , we think the backscattering effect 
will only change the prefactors of the resonance peak at ui — 002 = \Qvtl\- 

In conclusion, we have studied the RS from a TL liquid. If we treat the electron-photon 
matrix element as constant i.e. independent of wave- vector, the RS efficiency can not be 
simplified to be proportional to the structure factor S{qi — q2,^^i — ^^2)- Generally it has 
some complicated structure which is sensitive both to the coupling between the electron and 
core hole and the electron-electron interaction strength of the TL liquid. If we take into 
account the ^-dependence of M^, in the resonance region, the RS efficiency does have some 
part of the contribution which is proportional to the structure factor S{qi — q2, uJi — UJ2). We 
found the total secondary emission can not be rigorously separated into Hot Luminescence 
and Raman Scattering. 
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FIGURE CAPTIONS 



Fig.l Total HL transition rate Wp'^ + as a function of Ui (incident photon fre- 

quency) and UJ2 (scattered photon frequency) at g+ = 0.8, g_ = 0.4. The unit for the x, y 
axes is \Q\vvf. The positions of the two peaks are at CU1+CU2 — ±Qvtl- The singular line is 
at LOi — u;2 — \Q\vtl- Inset (a): The cross section at loi — u;2 — \Q\vtl + 0"^; Inset (b): The 
cross section at Ui = ^/2\Q\vtl- 

Fig. 2 The divergent region and convergent region of Wp. 
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